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Abstract: We have constructed models for a series of platinum-DNA adducts that represent 
the binding of two agents, [Pt 2 (DTBPA)Cl 2 ](H) and [Pt 2 (TPXA)Cl 2 ](II), to DNA via 
inter- and intra-strand cross-linking, and carried out molecular dynamics simulations and 
DNA conformational dynamics calculations. The effects of trans- and c/s-configurations of 
the centers of these di-nuclear platinum agents, and of different bridging linkers, have been 
investigated on the conformational distortions of platinum-DNA adducts formed via 
inter- and intra-strand cross-links. The results demonstrate that the DNA conformational 
distortions for the various platinum-DNA adducts with differing cross-linking modes are 
greatly influenced by the difference between the platinum-platinum distance for the 
platinum agent and the platinum-bound N7-N7 distance for the DNA molecule, and by the 
flexibility of the bridging linkers in the platinum agent. However, the effects of 
trans /cis-configamtions of the platinum-centers on the DNA conformational distortions in 
the platinum-DNA adducts depend on the inter- and intra-strand cross-linking modes. In 
addition, we discuss the relevance of DNA base motions, including opening, shift and roll, to 
the changes in the parameters of the DNA major and minor grooves caused by binding of the 
platinum agent. 
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1. Introduction 

Platinum-based drugs, such as cisplatin, carboplatin and oxaliplatin, are agents that have been used 
effectively to treat various types of human cancer, including genitourinary, colorectal and non-small 
cell lung cancers [1-5]. Various investigations have revealed that DNA is most likely the biological 
target associated with the antitumor activity of these drugs [6,7]. 

The platinum at the center of these agents binds to nucleophilic sites of DNA (the N7 atoms of 
guanosine [G] and adenosine [A] bases), and these interactions result in distortions of the DNA 
conformation that are critical to the antitumor effects [8-10]. However, the use of these platinum-based 
agents has often been limited by the development of tumor resistance [11-14]. Therefore, considerable 
effort has been made to develop new generations of platinum agents that do not share these 
limitations [15-17]. Polynuclear platinum agents represent a novel class of antitumor metallodrugs, in 
which two or three platinum centers are bridged by polyamines of variable length. These agents bind to 
DNA by forming various long-range intra-strand and inter-strand cross-links [18-20]. This unique 
DNA-binding feature is believed to be responsible for distortion of the DNA conformation, enhanced 
cytotoxicity, and effective prevention of resistance to the antitumor drug [21-23]. 

The di-nuclear platinum(II) agents, [Pt 2 (DTBPA)Cl 2 ](II) (where DTBPA = 2,20-(4,ll-dimethyl- 
1,4,8,1 l-tetraazacyclotetradecane-l,8-diyl)bis (A^-(2-(pyridin-2-yl)ethyl)acetamide)) and [Pt 2 (TPXA)Cl 2 ](II) 
(where TPXA = A^,A^,A^ , ,A^'-tetra(2-pyridylmethyl)-m-xylylene diamine), are two examples of 
polynuclear platinum drugs that are anticancer therapeutic agents [24,25]. In the [Pt 2 (DTBPA)Cl 2 ](II) 
agent, two platinum centers are bridged by a cyclam ring, and each platinum center is coordinated by 
one chloride and three nitrogens in, respectively, the cyclam ring, acetamide and pyridyl groups 
(see Chart 1) [24]. In the [Pt 2 (TPXA)Cl 2 ](II) agent, two platinum centers are bridged by a 
benzene linker, and each platinum center is coordinated by one chloride and three nitrogens in the 
bis(2-pyridylmethyl)amine (BPA) moiety (see Chart 1) [25]. Previous experiments have demonstrated 
that both agents exhibit significantly enhanced cytotoxicity against cancers compared with 
cisplatin [19-22]. Whereas cisplatin binds to DNA by forming short-range cross-links [26-28], the two 
di-nuclear platinum(II) agents, which have a cyclam ring or benzene group as a bridged linker, can 
bind to DNA by forming long-range inter-/intra- strand cross-links [24,25]. Moreover, the features of 
the multi-ring ligands in these two di-nuclear platinum(II) agents may facilitate platinum-DNA 
interactions through pre-associative electrostatic and H-bonding interactions [29-32]. Taking into 
account the platinum-platinum distances of 9.1-9.7 A in the [Pt 2 (DTBPA)Cl 2 ](II) and 
[Pt 2 (TPXA)Cl 2 ](II) agents, and the distance of about 3.4 A between two neighboring base pairs in 
B-DNA [24,25,33,34], it would be predicted that in the platinum-DNA adducts, the two platinum 
centers of [Pt 2 (DTBPA)Cl 2 ](II) or [Pt 2 (TPXA)Cl 2 ](II) would be capable of spanning 3-5 base pairs. 
Theoretically, Mori et al. 9 using ab initio fragment molecular orbital-based molecular dynamics 
(FMO-MD) simulations, presented that the Pt-Cl bond fluctuation and charge distribution for cis- and 
^rans-platinum are important for activating their anticancer functions [35,36]. However, few theoretical 
studies devoted to the properties of the [Pt 2 (DTBPA)Cl 2 ](II) and [Pt 2 (TPXA)Cl 2 ](II) agents binding to 
DNA molecules. 



Int. J. Mol Sci. 2013, 14 



19558 



Chart 1. X-ray structures for (a) [Pt 2 (DTBPA)Cl 2 ](II); and (b) [Pt 2 (TPXA)Cl 2 ](II); 
(c) the 18-mer duplex DNA sequence and the six platinum-DNA cross-linking modes for 
[Pt 2 (DTBPA)Cl 2 ](II) studied in this work, along with their assigned names. Three 
types of intra-strand cross-links (I: 26.28DTBPA-DNA; II: 4.7DTBPA-DNA; and 
III: 3.7DTBPA-DNA), and three types of inter-strand cross-links (IV: 7.28DTBPA-DNA; 
V: 6.28DTBPA-DNA; and VI: 7.26DTBPA-DNA) were studied. Similar nomenclature was 
used for [Pt 2 (TPXA)Cl 2 ](II). 




(b) (c) 



The [Pt 2 (DTBPA)Cl 2 ](II) agent, with two ^rans -platinum centers and a large 14-membered ring 
linker (cyclam), and the [Pt 2 (TPXA)Cl 2 ](II) agent, with two c/s-platinum centers and a small 
6-membered ring linker (benzene), exhibit distinct structural features that may influence their DNA 
binding modes and DNA distortion properties. However, how the structures of these two platinum 
agents affect their interaction and binding with DNA is currently not fully understood at the molecular 
level. In the present work, we have used molecular dynamics simulations and DNA conformational 
dynamics analyses to study variations in the conformations of the [Pt 2 (DTBPA)Cl 2 ](II)-DNA and 
[Pt 2 (TPXA)Cl 2 ](II)-DNA adducts with different cross-linking modes. Six different DNA-adduct 
models for each platinum agent, with inter-/intra-strand cross-linking modes of 1-3, 1-4 and 1-5 base 
intervals, were investigated to examine the properties of the DNA distortions induced by the two 
platinum agents. 

2. Results 

The 12 platinum-DNA adduct models (six [Pt 2 (DTBPA)Cl 2 ](II)+DNA adducts: 
26.28DTBPA-DNA, 4 . 7DTBP A-DN A, 3 . 7DTBP A-DN A, 7.28DTBPA-DNA, 6.28DTBPA-DNA and 
7.26DTBPA-DNA; and six [Pt 2 (TPXA)Cl 2 ](II) + DNA adducts: 26.28TPXA-DNA, 4.7TPXA-DNA, 
3.7TPXA-DNA, 7.28TPXA-DNA, 6.28TPXA-DNA and 7.26TPXA-DNA) were studied by 
performing 50-ns unrestrained molecular dynamics simulations. The root-mean-square deviation 
(RMSD) values for all the backbone atoms, referenced to the corresponding starting structures, over 
12 trajectories for the various intra-/inter-strand cross-linking models of the platinum-DNA adduct, were 
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examined to determine if each system had attained equilibrium. A small RMSD value for a simulation 
is often considered to indicate that the system is in a stable state, whereas large RMSD values suggest 
the occurrence of large conformational changes. Plots of the RMSDs of the 12 platinum-DNA system 
simulations over time are shown in Figure 1. The calculated conformational free energies of the 
platinum-DNA systems as a function of time, during the last 10 ns of the simulation, were 
analyzed [5,37]. The corresponding results for the 4 . 7DTBPX A-DN A, 6.28DTBPA-DNA, 
4.7TPXA-DNA and 6.28TPXA-DNA models are shown in Figure Al (a, b, c and d, respectively). It may 
be seen from Figures 1 and Al that the platinum-DNA adduct systems reached equilibrium after 10 ns, and 
that their energies were found to be stable during the remainder of each simulation. Therefore, the 
trajectory analysis for the 12 systems extracted the equilibrated conformations between 10 ns and 50 ns 
of the simulation time, recording 20,000 snapshots at 2-ps time intervals for each trajectory. The 
analyses of the hydrogen bonds of the DNA base pairs and the DNA-groove parameters around the 
binding sites of the DNA were performed using the time-averaged structures for these models, and are 
shown in Table 1, Figures 2 and 3. 

Figure 1. Root-mean- square deviation (RMSD) values of all the backbone atoms with 
respect to the corresponding starting structures for: (a) simulations of B-DNA (black), 
26.28DTBPA-DNA (blue), 4.7DTBPA-DNA (red) and 3.7DTBPA-DNA (green); 
(b) simulations of B-DNA (black), 7.28DTBPA-DNA (blue), 6.28DTBPA-DNA (red) and 
7.26DTBPA-DNA (green); (c) simulations of B-DNA (black), 26.28TPXA-DNA (blue), 
4.7TPXA-DNA (red) and 3.7TPXA-DNA (green); (d) simulations of B-DNA (black), 
7.28TPXA-DNA (blue), 6.28TPXA-DNA (red) and 7.26TPXA-DNA (green). 
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Figure 2. Hydrogen bond occupancies of the base pairs and the base-pair steps for the 
various platinum-DNA adducts (shown as circles for [Pt2(DTBPA)Cl2(II)] and squares for 
[Pt 2 (TPXA)Cl 2 (II)]). (a) 26.28DTBPA-DNA (green), 4.7DTBPA-DNA (red) and 
3.7DTBPA-DNA (violet); (b) 7.28DTBPA-DNA (cyan), 6.28DTBPA-DNA (orange) and 
7.26DTBPA-DNA (purple); (c) 26.28TPXA-DNA (magenta), 4.7TPXA-DNA (black) and 
3.7TPXA-DNA (blue); (d) 7.28TPXA-DNA (navy), 6.28TPXA-DNA (pink) and 
7.26TPXA-DNA (olive). 
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Figure 3. (a) The width and depth of the major groove of the time-averaged DNA 
conformation are shown as blue line (1) and green line (2), respectively, and the width and 
depth of the minor groove are shown as blue line (3) and green line (4), respectively; (b) The 
bend angles of the time-averaged DNA conformation. 
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Table 1. The widths and depths of the minor and major grooves and the bend angles for the 
time-averaged structures of the DNA conformations for the Pt 2 (DTBPA)/Pt 2 (TPXA)-DNA 
adducts and B-DNA. 



Model 


Major 
groove 
width (A) 


Major 
groove 
depth (A) 


Minor 
groove 
width (A) 


Minor 
groove 
depth (A) 


Bend (°) 


26.28DTBPA-DNA 


12.47 


3.28 


5.65 


4.93 


35.48 


4.7DTBPA-DNA 


13.72 


3.22 


5.95 


5.03 


30.70 


3.7DTBPA-DNA 


13.09 


4.94 


6.76 


4.15 


48.62 


7.28DTBPA-DNA 


14.85 


2.59 


6.21 


5.08 


53.81 


6 . 2 8DTBP A-DN A 


15.76 


2.46 


6.12 


4.29 


49.91 


7.26DTBPA-DNA 


13.66 


4.62 


5.03 


4.85 


61.94 


26.28TPXA-DNA 


12.27 


5.70 


5.91 


5.04 


27.20 


4.7TPXA-DNA 


12.43 


7.01 


6.81 


3.71 


26.18 


3.7TPXA-DNA 


10.52 


6.18 


6.53 


3.27 


49.15 


7.28TPXA-DNA 


12.95 


6.89 


5.52 


3.83 


30.58 


6.28TPXA-DNA 


11.98 


5.73 


6.04 


4.28 


27.15 


7.26TPXA-DNA 


9.73 


6.40 


6.35 


3.97 


50.57 


B-DNA 


11.81 


5.84 


5.87 


4.64 


15.02 



2.1. Conformational Characteristics and Flexibility of [Pt2(DTBPA)Ch] (II) and [Pt2(TPXA)Ch] (II) 

The [Pt2(DTBPA)Cl2](II) agent consists of a large cyclam ring linker and two ^rans-platinum(II) 
centers, each of which is coordinated by one chloride and one nitrogen each from the cyclam ring 
acetamide and pyridyl groups. The [Pt2(TPXA)Cl2](II) agent consists of a small benzene ring linker 
and two c/s-platinum(II) centers, each of which is coordinated in a square-planar geometry by one 
chloride and three nitrogen atoms from each bis(2-pyridylmethyl)amine (BPA) moiety. The 
conformational characteristics of the DNA binding processes for these two platinum(II) antitumor 
agents are expected to differ, due to the different ligands and linkers of the platinum(II) centers [24,25]. 
In order to compare the conformational flexibility of [Pt 2 (DTBPA)Cl 2 ](II) and [Pt 2 (TPXA)Cl 2 ](II), the 
average platinum-platinum distances of these platinum-DNA adducts were evaluated before and after 
simulation, using the conformation overlap method. Figure 4 shows the full-overlap conformations of 
[Pt 2 (DTBPA)Cl 2 ](II) and [Pt 2 (TPXA)Cl 2 ](II) between the crystal structures and averaged structures 
during the 50-ns simulations. It may be seen from Figure 4 that the average platinum-platinum distance 
values for [Pt 2 (DTBPA)Cl 2 ](II) and [Pt 2 (TPXA)Cl 2 ](II) change, during simulation, from 9.70 A and 
9.07 A in the crystal structures to 9.63 A and 6.98 A in the platinum-DNA adducts. Such small 
rangeability of the platinum-platinum distance for [Pt2(DTBPA)Cl2](II), which has a large-sized linker, 
during the simulations demonstrates that it has a more rigid conformation than [Pt2(TPXA)Cl2](II), 
which has a smaller- sized linker. 
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Figure 4. Overlap conformations between the crystal structures (white) and the averaged 
simulated structures of (a) [Pt 2 (DTBPA)Cl 2 ](II) (gray); and (b) [Pt 2 (TPXA)Cl 2 ](II) (gray). 
The platinum atoms are shown as white in the crystal structures, and light cyan in the 
averaged structures. 




(b) 



2.2. Relationship between the Platinum-Platinum Distances in the Platinum(II) Agents and the N7-N7 
Distances in the Pt 2 (DTBPA)-DNA and Pt 2 (TPXA)-DNA Adducts 

To address the effects of the platinum-platinum distance and platinum-bound N7-N7 distance on the 
extent of the distortion of the DNA conformation in the platinum-DNA adduct, the overlaps of the 
averaged structures for intra- 1,3-GG cross-linking 26.28DTBPA-DNA, intra- 1,4-GG cross-linking 
4.7DTBPA-DNA and intra- 1,5-GG cross-linking 3.7DTBPA-DNA were compared with a bare DNA 
molecule (Figure 5). The platinum-bound N7-N7 distance in the three models was used as a test 
parameter for the degree of the DNA conformational distortion. The platinum-bound N7-N7 distances 
in DNA and Pt-Pt distances in the platinum agents for the three models and for the B-DNA molecule are 
shown in Table 2. Of the three models, the platinum-platinum distance in the [Pt 2 (DTBPA)Cl 2 ](II) 
agent most closely corresponded to the platinum-bound N7-N7 distance in the intra- 1,4-GG 
cross-linking 4.7DTBPA-DNA. It may be seen from Table 2 that the difference in the N7-N7 distance 
between the B-DNA molecule and the averaged platinum-bound adduct was 2.22 A for the 
intra-l,3-GG cross-linking 26.28DTBPA-DNA, and 3.76 A for the intra-l,5-GG cross-linking 
3.7DTBPA-DNA adducts, but only 0.38 A for the intra-l,4-GG cross-linking 4.7DTBPA-DNA. 
Moreover, as shown in Figure 2a, the sum of the percentages of the hydrogen bonds destroyed at the 
bases, in comparison with the B-DNA molecule, was 262.55% for 26.28DTBPA-DNA (A10, 89.79%; 
Cll, 95.13%; and A12, 77.63%) and 404.57% for 3.7DTBPA-DNA (A2, 37.87%; A3, 100%; 
G4, 92.66%; A5, 75%; and A6, 99.04%), but only 202.05% for 4.7DTBPA-DNA (G4, 69.88%; 
A5, 57.79%; and A6, 74.38%). The differences in the bend angles of the DNA molecule toward the 
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minor groove between the B-DNA molecule and each adduct were 20.46° for 26.28DTBPA-DNA, 
33.60° for 3.7DTBPA-DNA, and only 15.68° for 4.7DTBPA-DNA (see Table 1 and Figure 3). As would 
be expected, Figure 5 indicates visually that the distortion of the average DNA conformation, compared 
with the B-DNA molecule, was also larger for both 26.28DTBPA-DNA and 3.7DTBPA-DNA models 
than for 4.7DTBPA-DNA model. Therefore, large differences between the platinum-platinum distance 
in the platinum agent and the platinum-bound N7-N7 distance in the DNA molecule can cause 
substantial distortions of the DNA conformation during the process of platinum-DNA adduct formation. 
Similar results were observed for the [Pt2(TPXA)Cl2](II) agent and the inter-strand cross-linking 
models, except for the fact that the effect on the conformationally-rigid [Pt2(DTBPA)Cl2](II) agent was 
more obvious than that on [Pt2(TPXA)Ci2](II), which has a more flexible conformation. 

Figure 5. The averaged centroid structures of the platinum-DNA adduct and an undamaged 
B-DNA for: (a) 26.28DTBPA-DNA; (b) 4.7DTBPA-DNA; and (c) 3.7DTBPA-DNA. 
[Pt 2 (DTBPA)Cl 2 ](II) is shown as pink, bound-DNA as blue, and B-DNA as gray. 




(a) (b) (c) 

Table 2. Differences between the platinum-binding atoms of the base pairs (N7-N7) of bare 
DNA and the Pt 2 (DTBPA)-DNA adduct, before and after simulation (A(Pt-N7)); and 
between the Pt-Pt distance in the [Pt 2 (DTBPA)Cl 2 ](II) agent and N7-N7 distance in bare 
DNA (A(N7-N7)). 





Pt-Pt 


N7-N7 distance 






Model 


distance 

(A) 


B-DNA (A) 


Platinum-DNA 
adduct (A) 


A(Pt-N7) (A) 


A(N7-N7) (A) 


3.7DTBPA-DNA 


9.70 


14.81 


11.05 


5.11 


3.76 


4.7DTBPA-DNA 


9.70 


11.60 


11.22 


1.90 


0.38 


26.28DTBPA-DNA 


9.70 


7.63 


9.85 


2.07 


2.22 
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2.3. Relationship between the Flexibilities of the fPt 2 (DTBPA)Cl 2 J(II) and [Pt 2 (TPXA)Cl 2 ](II) Agents 
and Their DNA Binding Interactions 

Based on the fact that the [Pt 2 (DTBPA)Cl 2 ](II) agent has more rigidity than the [Pt 2 (TPXA)Cl 2 ](II) 
agent, a comparison was made of the DNA conformational distortions caused by the rigid 
[Pt 2 (DTBPA)Cl 2 ](II) agent and the flexible [Pt 2 (TPXA)Cl 2 ](II) agent, using the 4.7DTBPA-DNA and 
4.7TPXA-DNA adducts. Figure 6a shows the overlaps of the averaged structures for these two models, 
which indicate that the distortion of the DNA conformation for the 4.7DTBPA-DNA model was larger 
than that for the 4.7TPXA-DNA model. The sum of the percentages of the hydrogen bonds destroyed at 
bases G4, A5 and A6, in comparison with the B-DNA molecule, was 202.05% for the 4.7DTBPA-DNA 
adduct, and only 120.28% for the 4.7TPXA-DNA adduct (see Figure 2a,c). The average deviation of the 
major groove width, depth and bend angle were also much larger for the 4.7DTBPA-DNA adduct 
(16.17% [from 11.81 A to 13.72 A], 44.86% [from 5.84 A to 3.22 A] and 104.39% [from 15.02° to 
30.70°], respectively) than for the 4.7TPXA-DNA adduct (6.09% [from 11.81 A to 12.53 A], 14.21% 
[from 5.84 A to 5.01 A] and 74.30% [from 15.02° to 26.18°], respectively) (see Table 1 and Figure 3). 
Similar results for the other inter-/intra-strand cross-linking models were observed. Consequentially, the 
more rigid platinum adducts may lead to obvious distortion of the DNA backbone. 

Figure 6. Overlap of the averaged structures for B-DNA (gray), [Pt 2 (DTBPA)Cl 2 ](II) and 
[Pt 2 (TPXA)Cl 2 ](II) (pink), and: (a) 4.7DTBPA-DNA (blue) and 4.7TPXA-DNA (cyan); 
(b) 4.7DTBPA-DNA (blue) and 6.28DTBPA-DNA (royal); and (c) 4.7TPXA-DNA (cyan) 
and 6.28TPXA-DNA (green). 




(a) (b) (c) 
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2.4. Relationship between the trans- [Pt 2 (DTBPA)Cl 2 ] 1 (II)/c\s-[Pt 2 (TPXA)Cl 2 ] } (II) Configurations and 
the DNA Inter-/Intra-Strand Cross-linking Modes 

To investigate the effects of *ra^-[Pt 2 (DTBPA)Cl 2 ](II) and c/s-[Pt 2 (TPXA)Cl 2 ](II) conformations 
on the distortion of the DNA conformation in the intra-/inter-strand cross-linking modes of the 
platinum-DNA adducts, comparisons were made of the DNA conformational distortions between the 
B-DNA molecule and each of the intra- 1,4-GG cross-linking modes of the 4.7DTBPA-DNA and 
4.7TPXA-DNA adducts, and the inter- 1,4-GG cross-linking modes of the 6.28DTBPA-DNA and 
6.28TPXA-DNA adducts. Figures 6b,c show the overlaps of the averaged structures for these four 
adducts, indicating that for the trans-type platinum agent, the distortion of the DNA conformation (in 
comparison with the B-DNA molecule) was larger for the inter-strand cross-linking 6.28DTBPA-DNA 
model than for the intra-strand cross-linking 4.7DTBPA-DNA model; in contrast, for the cis-type 
platinum agent, the DNA conformational distortion was larger for the intra-strand cross-linking 
4.7TPXA-DNA model than for the inter-strand cross-linking 6.28TPXA-DNA model. The sum of the 
hydrogen bonds destroyed at the bases for / L ran l s"-[Pt 2 (DTBPA)Cl 2 ](II) was 282.06% for the inter-strand 
cross-linking 6.28DTBPA-DNA adduct (A5, 10.91%; A6, 58.14%; G7, 48.16%; T8, 81.05%; 
C9, 78.43%; and A10, 5.37%), much larger than the value of 202.05% for the intra-strand cross-linking 
4.7DTBPA-DNA adduct (at the G4, A5 and A6 bases) (see Figure 2a,b). In contrast, for the 
c/s-[Pt 2 (TPXA)Cl 2 ](II) agent, the sum of the hydrogen bonds destroyed at bases G4, A5 and A6 was 
120.28%) for the intra-strand cross-linking 4.7TPXA-DNA adduct, much larger than the value of 55.04%> 
at only one base (A10) for the inter-strand cross-linking 6.28TPXA-DNA adduct (see Figure 2c,d). 
Similarly, for the £rans-[Pt 2 (DTBPA)Cl 2 ](II) agent, the average deviations of the major groove width, 
depth and bend angle, in comparison with the B-DNA molecule, were larger for the inter-strand 
cross-linking 6.28DTBPA-DNA adduct (33.44% [from 11.81 A to 15.76 A], 57.87% [from 5.84 A to 
2.46 A] and 232.29%) [from 15.02° to 49.91°], respectively) than for the intra-strand cross-linking 
4.7DTBPA-DNA adduct (16.17% [from 11.81 A to 13.72 A], 44.86% [from 5.84 A to 3.22 A] and 
104.39% [from 15.02° to 30.70°], respectively). In contrast, for the c/s-[Pt 2 (TPXA)Cl 2 ](II) agent, the 
corresponding values were larger for the intra-strand cross-linking 4.7TPXA-DNA adduct (6.09%) [from 
1 1.81 A to 12.53 A], 14.21% [from 5.84 A to 5.01 A] and 74.30% [from 15.02° to 26.18°], respectively) 
than for the inter-strand cross-linking 6.28TPXA-DNA adduct (1.44% [from 11.81 A to 11.98 A], 
1.88% [from 5.84 A to 5.73 A] and 80.76% [from 15.02° to 27.15°], respectively) (see Table 1 and 
Figure 3). Similar results were observed for other inter-/intra-strand cross-linking models of these two 
platinum agents. Taken together, these data suggest that the trans-type platinum agents facilitate greater 
DNA conformational distortion through an inter-strand rather than an intra-strand cross-linking 
mode, whereas the cis-type platinum agents create greater DNA conformational distortion through an 
intra-strand rather than an inter- strand cross-linking mode [38-40]. Furthermore, platinum agents 
with trans-type DNA cross-links have been shown to have high antitumor activity in 
cytotoxicity experiments [41]. 
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2.5. Configuration Characteristics of Two Additional Platinum Agents, tmns-[Pt2(DPZM)Cl2] (II) and 
cis-[Pt 2 (DPZMCH2)Cl 2 ](II), Binding to DNA 

To confirm generalizability of these relationships in [Pt 2 (DTBPA)Cl 2 ](II) and [Pt 2 (TPXA)Cl 2 ](II) 
agents, the two additional dinuclear platinum agents, £rans-[Pt 2 (DPZM)Cl 2 ](II) (where 
DPZM = {(NH 3 )} 2 u-4,4'-dipyrazolylmethane) and c^-[Pt 2 (DPZMCH2)Cl 2 ](II) (where 
DPZMCH2 = {(NH 3 )} 2 u-CH 2 CH 2 -4,4'-dipyrazolylmethane) binding to the same DNA molecule were 
simulated to study the characteristics of DNA distortion [42] (see Chart Al). In the [Pt 2 (DPZM)Cl 2 ](II) 
agent, two trans -platinum centers are bridged by one acyclic alkane linker (-CH 2 -); each platinum 
center is coordinated by one chloride and three nitrogens in, respectively, ammonia, ammonia and 
pyrazole groups. To improve the flexibility of the [Pt 2 (DPZM)Cl 2 ](II) agent, a relatively flexible linker 
(-CH 2 -CH 2 -CH 2 -) between two c/s-platinum centers is introduced to the [Pt 2 (DPZMCH2)Cl 2 ](II) 
agent. Six inter-/intra-strand cross-linking modes of 26.28DPZM-DNA, 4.7DPZM-DNA, 
3.7DPZM-DNA, 6.28DPZM-DNA, 4 . 7DPZMCH2-DN A and 6.28DPZMCH2-DNA were built and 
simulated by using the same methods described in the "models and methods" section. The 
corresponding results for these six additional models were shown in Tables Al and A2, Figures A2-A6. 
It can be seen from Tables Al and A2, Figures A3 and A5 that the distortion of the average DNA 
conformation, compared with the B-DNA molecule, was larger for both 26.28DPZM-DNA and 
3.7DPZM-DNA models than for the 4.7DPZM-DNA model due to the differences between the Pt-Pt 
distance of platinum agent and the N7-N7 distances of the B-DNA molecule. For example, the sum of 
the percentages of the hydrogen bonds destroyed was 400.06% for the 3.7DPZM-DNA model with this 
distance difference of 5.06 A versus 319.72% for the 4.7DPZM-DNA model with that of 1.85 A. The 
results further confirm that large differences between the platinum-platinum distance in the platinum 
agent and the platinum-bound N7-N7 distance in the DNA molecule can cause substantial distortions of 
the DNA conformation discussed above for the two [Pt 2 (DTBPA)Cl 2 ](II) and [Pt 2 (TPXA)Cl 2 ](II) 
agents. It can be seen from Table Al, Figures A3, A4 and A6 that the distortion of the DNA 
conformation for the rigid 4.7DPZM-DNA model with the deviation of DNA major groove width and 
depth of 9.65% and 30.06% was larger than that for the flexible 4.7DPZMCH2-DNA model with those 
of 5.08%) and 13.35%, which confirms that the rigid platinum adduct may lead to more obvious 
distortion of the DNA backbone than the flexible platinum adduct. Moreover, for the 
/ L ran l s"-[Pt 2 (DPZM)Cl 2 ](II) agent, the distortion of the DNA conformation (in comparison with the 
B-DNA molecule) was larger for the inter-strand cross-linking 6.28DPZM-DNA model than for the 
intra-strand cross-linking 4.7DPZM-DNA model; in contrast, for the c/s-[Pt 2 (DPZMCH2)Cl 2 ](II) agent, 
the DNA conformational distortion was larger for the intra-strand cross-linking 4.7DPZMCH2-DNA 
model than for the inter-strand cross-linking 6.28 DPZMCH2-DNA model, which further verifies that 
the trans-type platinum agents facilitate greater DNA conformational distortion through an inter-strand 
rather than an intra-strand cross-linking mode (see Table Al, Figures A3 and A6). 
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3. Discussions 

3.1. Principal Component Analysis of the Major Conformational Dynamics 

Principal component analysis was used to investigate the trajectories from the corresponding 
simulations to examine the dominant dynamic motions of the DNA in the 
12 Pt 2 (DTBPA)/Pt 2 (TPXA)-DNA adducts studied. The first three principal components (PCI, PC2 and 
PC3) described about 90% of the essential modes of the dynamics for these adducts (Table 3). It was 
found that the first three components of the conformational motions roughly corresponded to a 
superposition of bending, unwinding and twisting motions. Visual analyses of averaged DNA structures 
from the trajectories support the dominant motions identified by PC A. The overall centroid structures 
for these adducts have been analyzed, and are shown in Figure A7. It may be seen that compared with 
the undamaged B-DNA conformation, there were obvious differences between the various 
cross-linking models in the DNA conformations of the platinum-DNA adducts (including bending, 
unwinding and twisting motions). In particular, the PCI occupancy time percentages for binding of 
[Pt 2 (DTBPA)Cl 2 ](II) to DNA were larger than those for [Pt 2 (TPXA)Cl 2 ](II), which would be expected 
to lead to larger DNA bending angles for the rigid [Pt 2 (DTBPA)Cl 2 ](II) agent than for the flexible 
[Pt 2 (TPXA)Cl 2 ](II) agent (with respect to the undamaged B-DNA molecule) (see Table 3). The 
unwinding motions, in the vicinity of the platinum sites, for the Pt 2 (DTBPA)/Pt 2 (TPXA)-DNA adducts 
also efficiently induce DNA groove conformational changes; the first three motions take place 
simultaneously and are irreversible. 



Table 3. Occupancy time percentages for the first three principal components during 
the simulations for each of the studied Pt 2 (DTBPA)/Pt 2 (TPXA)-DNA adducts (PCsum 
represents the sum of PCI, PC2 and PC3). 



Model 


PCI 


PC2 


PC3 


PCsum 


26.28DTBPA-DNA 


77.60 


7.74 


7.07 


92.41 


4.7DTBPA-DNA 


66.56 


22.98 


1.98 


91.52 


3.7DTBPA-DNA 


48.93 


32.79 


12.50 


94.22 


7.28DTBPA-DNA 


67.36 


20.31 


5.89 


93.56 


6.28DTBPA-DNA 


67.37 


22.45 


3.12 


92.94 


7.26DTBPA-DNA 


54.53 


21.73 


16.58 


92.84 


26.28TPXA-DNA 


60.40 


26.49 


4.22 


91.11 


4.7TPXA-DNA 


58.00 


24.51 


7.66 


90.17 


3.7TPXA-DNA 


63.53 


19.50 


8.80 


91.83 


7.28TPXA-DNA 


62.46 


19.07 


8.45 


89.98 


6.28TPXA-DNA 


40.14 


331.45 


116.72 


888.31 


7.26TPXA-DNA 


50.27 


24.79 


15.20 


90.26 



3.2. DNA Conformational Dynamics in the Platinum-DNA Adducts 

The [Pt 2 (DTBPA)Cl 2 ](II) and [Pt 2 (TPXA)Cl 2 ](II) agents disturb the DNA conformation through the 
covalent binding of the platinum centers to the N7 atoms of the DNA bases, and simultaneously, the 
locating of the corresponding large ligand of the agent to the major groove. The processes of DNA 
conformational dynamics result in changes in the DNA groove parameters. The frequency distributions 
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of the DNA helical dynamics parameters have been analyzed, and are shown in Figure 7 for the 
4.7DTBPA-DNA adduct, and in Figure A8 for the other 1 1 models of Pt 2 (DTBPA)/Pt 2 (TPXA)-DNA. In 
general, groove width and depth were measured mainly from the opening, twist and rise motions, and the 
shift motion, respectively, of the DNA base pair helixes or steps, as defined by CURVES 
analysis [43-45]. For example, for the 4.7DTBPA-DNA model (see Figure 7), the opening parameters 
of the DNA base pair helixes at connection sites G4:C33, A5:T32 and A6:T31 were changed — 90°, 
— 80° and — 80°, respectively, away from the B-DNA molecule (compared with -0° for the B-DNA 
molecule); the deviations in the twist angle values were -20°, -30° and — 35° for the G4- A5, A5 • A6 and 
A6G7 base pair steps, with respect to the value of -35° for the B-DNA molecule; and the rise value 
deviations were -1.5 A, — 1 A and ~5 A for the G4- A5, A5- A6 and A6 G7 base pair steps, respectively, 
with respect to the value of -3.5 A for the B-DNA molecule. This results in a widening of the DNA 
major groove by 1.91 A around base pairs G4:C33-G7:C30 (see Table 1 and Figure 8a), and the 
destruction of hydrogen bonds at base pairs G4:C33, A5:T32 and A6:T31 (see Figure 2a). The 
deviations in the shift values were ~8 A, ~2 A and — 1 A for the G4A5, A5A6 and A6G7 base pair 
steps, respectively, compared with ~0 A for the B-DNA molecule (see Figure 7); this results in a shoaling 
of the DNA major groove by 2.62 A around base pairs G4:C33-G7:C30 (see Table 1 and Figure 8b). In 
addition, the DNA bend was measured by the twist, tilt and roll motions for DNA base pair steps. For 
example, for the 4.7DTBPA-DNA model (see Figure 7), the deviations in the values of the roll angles 
were — 40°, — 25°, and -40° for the G4-A5, A5A6 and A6G7 base pair steps, with respect to -0° for 
the B-DNA molecule; this contributes to a significant bending of the DNA helix toward the minor 
groove of 15.68°. Similar results were found for the other 1 1 models of Pt 2 (DTBPA)/Pt 2 (TPXA)-DNA 
(see Figure A9). Taken together, these data indicate that binding of the two types of platinum agent to the 
DNA molecule greatly affects the DNA groove parameters via dynamic changes in DNA conformation 
around the connected base pairs, leading to difficulties in DNA-repair as well as gene transcription. 

Figure 7. Selected frequency distributions of the representative DNA duplex base-pair 
helical/step parameters for the central binding base pairs of the 4.7DTBPA-DNA model 
(red line), compared with B-DNA (black line). 
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Figure 8. Major groove width (a) and depth (b) for the time-averaged structure of the DNA 
conformation of the 4.7DTBPA-DNA model (red line with circles), compared with 
B-DNA (black line with squares). 



20 



16 



i 



o) 12 
o 



-m- B-DNA 
• 4.7DTBPA-DNA 




10 



7 8 9 10 11 12 13 14 15 16 
Base-Pair number 



(a) 



B-DNA 
4.7DTBPA-DNA 




7 8 9 10 11 12 13 14 15 16 
Base-Pair number 



(b) 



Int. J. Mol. Sci. 2013, 14 



19570 



4. Models and Methods 

4.1. Starting Structures 

The starting structures of all the platinum-DNA adducts contained an 18-mer duplex DNA sequence, 
5 f -d(GAAGAAGTCACAAAATGT)-3 f , that was chosen on the basis of the breast/ovarian 
cancer-susceptibility gene, BRCA1 (Genbank #U14680) [10], and the [Pt 2 (DTBPA)Cl 2 ](II) or 
[Pt2(TPXA)Cl2](II) agent selected from X-ray data [24,25]; see Chart 1. Previous studies have 
suggested that the two platinum center atoms may bind to the N7 atoms of the G/A bases of DNA, and 
might span 3-5 base pairs with intra-/inter-strand cross-links. Therefore, three types of intra-strand 
cross-link were constructed and simulated for [Pt2(DTBPA)Cl2](II), with 1232 atoms: an intra-l,3-GG 
cross-linking mode (designated 26.28DTBPA-DNA), an intra- 1,4-GG cross-linking mode (designated 
4.7DTBPA-DNA), and an intra- 1,5- AG cross-linking mode (designated 3 . 7DTBP A-DN A) ; and three 
types of inter-strand cross-link were constructed and simulated for [Pt2(DTBPA)Cl2](II), with 
1232 atoms: an inter- 1,3-GG cross-linking mode (designated 7.28DTBPA-DNA), an inter- 1,4-AG 
cross-linking mode (designated 6.28DTBPA-DNA), and an inter- 1,5-GG cross-linking mode 
(designated 7.26DTBPA-DNA). The same protocol (with similar nomenclature) was applied to the 
[Pt 2 (TPXA)Cl 2 ](II) agent. Three of these 12 adduct models, 26.28TPXA-DNA, 4.7TPXA-DNA and 
7.28TPXA-DNA, have been previously examined using denaturing polyacrylamide gels and 
MALDI-TOF MS analysis [25]. The initial coordinates of the 12 platinum(II)-DNA adducts used in 
our simulations were generated by docking the platinum agent in the major groove of the DNA. For 
each binding adduct, several different orientations of the platinum compound, coordinating with the 
major groove of DNA, were chosen as the starting structures for our molecular dynamics simulations, 
thereby providing a good test of whether our molecular dynamics simulations were capable of driving 
significantly distinct starting structures to a non-distinguishable one when the simulations reached 
equilibrium. Bump-checking was turned on to ensure that no overlapping atoms were produced during 
the docking process. To compare the differences between the DNA adduct conformations and an 
undamaged DNA molecule, a bare-DNA (B-DNA) molecule simulation was also performed, in which 
an idealized B-DNA was used as a starting structure for the simulation. Given that each strand of DNA 
has numerous phosphate groups, 32 Na + counterions were added to each system to achieve 
electroneutrality. The systems were explicitly solvated using the TIP3P water potential inside a box 
large enough to ensure that the solvent shell extended to 10 A in all directions. 

4.2. Force Field Parameter Preparation 

The atoms for the studied platinum compound, excluding the platinum atoms themselves, were 
generated using the ANTECHAMBER module in the AMBER9 program [46]. The force field 
parameters and the RESP charges of the platinum compound were referenced from previous 
works [47-51]. Other force field parameters for the platinum compound were generated from the gaff 
force field in the AMBER9 program. 
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4.3. Molecular Dynamics Simulations and Trajectory Analysis 

All molecular dynamics simulations were carried out using the AMBER9 package [46] with a 
classical AMBER parm99 [52,53] together with the parmbscO refinement [48,54] and gaff [55] force 
field parameters. Details of the molecular dynamics protocols are given in the Appendix. 

4.4. Principal Component Analysis 

Principal component analysis can be used to segregate large-scale correlated motions from random 
thermal fluctuations, thereby probing the essential dynamics of the system. Details of the analysis 
method are available in the Appendix. 

4.5. DNA Groove Parameter Analyses of Trajectories 

To investigate the distortion of the DNA, the frequency distributions (fraction of the time spent in 
each conformation) of the trajectories of the simulations for the models and a canonical B-DNA were 
calculated using the CURVES program [56]. To account for distortion of the whole DNA backbone, the 
overall bend, tilt and roll angles of the DNA time-averaged structures for the studied models were 
calculated from the CURVES outputs using MadBend [57]. Details of the calculation method are 
available in the Appendix. 

5. Conclusions 

Molecular dynamics simulations and DNA dynamics analyses for a series of 
[Pt 2 (DTBPA)Cl 2 ](II)-DNA and [Pt 2 (TPXA)Cl 2 ](II)-DNA adducts, with inter-/intra- strand 
cross-linking modes and 1-3, 1-4, 1-5 base intervals, were carried out to examine the distortions in 
the DNA double-helical structure. For each series of platinum agents, it was found that a larger 
difference between the platinum-platinum distance of the platinum agent and the platinum-bound 
N7-N7 distance of the DNA molecule was associated with a greater DNA conformational distortion. For 
the same cross-linking mode, the rigid [Pt 2 (DTBPA)Cl 2 ](II) agent induced a larger DNA conformational 
distortion than the flexible [Pt 2 (TPXA)Cl 2 ](II) agent. The drafts -platinum configuration of 
[Pt 2 (DTBPA)Cl 2 ](II) may facilitate conformational distortion of the platinum-DNA adduct through 
inter-strand cross-linking modes, whereas the c/s-platinum configuration of [Pt 2 (TPXA)Cl 2 ](II) may 
favor intra-strand cross-linking modes; this is consistent with previous studies [50-52]. DNA dynamics 
analysis revealed that the DNA base motions (opening, shifting, rising, rolling and twisting away from 
the DNA grooves) caused by the binding of the platinum agent induced widening of the DNA major and 
minor grooves and destruction of the hydrogen bonds of the DNA base pairs. These observations are 
validated for the other two *rara-[Pt 2 (DPZM)Cl 2 ](II) and c/s-[Pt 2 (DPZMCH2)Cl 2 ](II) agents binding to 
DNA. Our calculated results provide useful insights, in atomic detail, into how DNA conformations 
are affected by platinum agents containing different ligands and linkers, and may aid the experimental 
design of novel types of platinum-based anticancer agent. 
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Appendix 

Molecular Dynamics Simulation Protocols 

All molecular dynamics (MD) simulations were carried out using the AMBER9 package [37,58] 
with a classical AMBER parm99 [43,44] together with the parmbscO refinement [39,45] and gaff [46] 
force field parameters. The protocol for all MD simulations was as follows: (1) All systems were 
submitted to energy minimization to remove unfavorable contacts. Four cycles of minimizations were 
performed with 2500 steps in each minimization, with harmonic restraints ranging from 
100 kcal mol -1 A 2 , 25 kcal mol -1 A -2 to 10 kcal mol 1 A 2 on the DNA and platinum compound 
positions. Finally, 5000 steps of unrestrained minimization were carried out before the heating process. 
The cutoff distance used for non-bonded interactions was 10 A. The SHAKE algorithm [59] was used 
to restrain the bonds containing hydrogen. (2) Each energy-minimized structure was heated over 120 ps 
from 0 to 300 K (with a temperature coupling of 0.2 ps) in a constant volume, while the positions of 
DNA and the platinum compound were restrained with a small value of 10 kcal mol 1 A 2 . A constant 
volume was maintained during the processes. (3) An unrestrained equilibration was carried out for 200 ps 
under constant pressure and constant temperature conditions, before a trajectory was generated for a 
further simulation. The temperature and pressure were allowed to fluctuate at around 300 K and 1 bar, 
respectively, with a corresponding coupling of 0.2 ps. For each simulation, an integration step of 2 fs 
was used. (4) Production runs of 50 ns were carried out using the same protocol. In each simulation, 
the time point at 200 ps after thermal equilibration was selected as the starting point for data collection. 
During the production run, 25,000 structures from each simulation were saved for post-processing by 
uniformly sampling the trajectory. (5) Finally, simulated annealing was performed for each system by 
allowing the system to cool from 300 K to 0 K during the course of a 1-ns simulation. The 
conformational free energies for all systems were calculated, using MM-PBSA energy analysis, from 
1000 snapshots extracted at 5-ps intervals during the last 10 ns of each simulation. The MM-PBSA free 
energy computations were carried out using the AMBER9 program. 

Principal Component Analysis 

Principal component analysis (PCA) was used to identify the dominant motions observed during a 
simulation by visual inspection; this was achieved by calculating the eigenvectors from the covariance 
matrix of a simulation and then filtering the trajectories along each of the different eigenvectors. As PCA 
reduces the dimensionality of the trajectory data, by using only the first several principal components we 
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can extract important dynamical features describing the collective and overall motions of the 
system [60-63]. The platinum-DNA adduct systems reached equilibrium after 10 ns; thus, the final 40 ns 
of the equilibrated MD simulations of the 20,000 trajectories for each of the platinum-DNA adducts was 
subjected to PCA analysis. The first three principal components described >90% of the essential modes 
of the dynamics for the platinum- adducted DNA complexes. The dynamics of the first three principal 
components were visualized using the Discovery Studio molecular graphics tool. 

DNA Groove Parameter Analyses of the Trajectories 

The PTRAJ module of the AMBER9 program was used to extract the produced conformations. 
These extracted snapshots were saved in the Protein Data Bank (PDB) format. Each nucleotide type 
was converted from the AMBER format to PDB format, and the resulting snapshots were submitted to 
the CURVES program [47]. The following CURVES parameters were extracted: global base-base 
parameters (shear, stretch, stagger, buckle, propeller and opening); and global inter-base-pair parameters 
(shift, slide, rise, tilt, roll and twist). Percentage occupancy distributions of the DNA helical parameters 
were calculated by normalizing the frequency distributions to 100%. 

The overall bend, tilt and roll angles of the DNA for the averaged structures of these adducts were 
calculated from the CURVES output using MadBend (Department of Chemistry and Courant Institute 
of Mathematical Sciences, New York University and Howard Hughes Medical Institute, New York, NY, 
USA http://www.biomath.nyu.edu/index/software/Madbend/index.html), developed by Strahs and 
Schlick [48]. This method evaluates the DNA curvature by summing the projected components of 
local base-pair step tilt and roll angles after adjusting the helical twist. Bends in the helical axis defined 
by a negative roll angle indicate bending toward the minor groove, while bends defined by a positive 
roll angle correspond to bending toward the major groove [48]. 



Table Al. The widths and depths of the minor and major grooves and the bend angles for 
the time-averaged structures of the DNA conformations for the various 
Pt 2 (DPZM)/Pt 2 (DPZMCH2)-DNA adducts and B-DNA. 



Model 


Major 
groove 
width (A) 


Major 
groove 
depth (A) 


Minor 
groove 
width (A) 


Minor 
groove 
depth (A) 


Bend 

(°) 


26.28DPZM-DNA 


13.27 


6.51 


4.67 


5.64 


43.91 


4.7DPZM-DNA 


12.95 


4.49 


6.42 


4.14 


30.88 


3.7DPZM-DNA 


13.11 


4.35 


6.23 


4.36 


34.94 


6.28DPZM-DNA 


15.40 


2.51 


4.16 


3.79 


48.73 


4.7DPZMCH2-DNA 


12.41 


6.62 


6.04 


3.98 


21.92 


6.28DPZMCH2-DNA 


11.47 


6.23 


6.51 


4.16 


19.24 


B-DNA 


11.81 


5.84 


5.87 


4.64 


15.02 
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Table A2. Differences between the platinum-binding atoms of the base pairs (N7-N7) of 
bare DNA and the Pt 2 (DPZM)-DNA adduct, before and after simulation (A(Pt-N7)); and 
between the Pt-Pt distance in the Pt2(DPZM)Cl2 agent and N7-N7 distance in bare 
DNA (A(N7-N7)). 





Pt-Pt 


N7-N7 distance 


A(Pt-N7) 

(A) 




Model 


distance 

(A) 


B-DNA (A) 


Platinum-DNA 
adduct (A) 


A(N7-N7) (A) 


3.7DPZM-DNA 


9.75 


14.81 


12.80 


5.06 


2.01 


4.7DPZM-DNA 


9.75 


11.60 


11.47 


1.85 


0.13 


26.28DPZM-DNA 


9.75 


7.63 


11.07 


2.12 


3.44 



Chart Al. X-ray structures for (a) ^ara-[Pt 2 (DPZM)Cl 2 ](II), and (b) c/s-[Pt 2 (DPZMCH2)Cl 2 ](II). 




(b) 

Figure Al. Conformational stability free energy (kcal-mol -1 ) and integrated distributions 
for (a) 4.7DTBPA-DNA (red); (b) 6.28DTBPA-DNA (blue), (c) 4.7TPXA-DNA (magenta) 
and (d) 6.28TPXA-DNA (green). 
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Figure Al. Cont. 
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Figure A2. RMSD values of all the backbone atoms with respect to the corresponding 
starting structures for: (a) simulations of B-DNA (black), 26.28DPZM-DNA (blue), 
4.7DPZM-DNA (red) and 3.7DPZM-DNA (green); (b) simulations of B-DNA (black) and 
6.28DPZM-DNA (red); (c) simulations of B-DNA (black) and 4 . 7DPZMCH2-DN A (red); 
(d) simulations of B-DNA (black) and 6.28DPZMCH2-DNA (red). 
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Figure A3. Hydrogen bond occupancies of the base pairs and the base-pair steps for the 
various platinum-DNA adducts (shown as circles for [Pt2(DPZM)Cl2(II)] and squares for 
[Pt 2 (DPZMCH2)Cl 2 (II)]). (a) 26.28DPZM-DNA (green), 4.7DPZM-DNA (red) and 
3.7DPZM-DNA (violet); (b) 6.28DPZM-DNA (orange); (c) 4 . 7DPZMCH2-DN A (black); 
(d) 6.28 DPZMCH2-DNA (pink). 
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Figure A3. Cont. 
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Figure A4. Overlap conformations between the crystal structures (white) and the averaged 
simulated structures of (a) [Pt 2 (DPZM)Cl 2 ](II) (gray); and (b) [Pt 2 (DPZMCH2)Cl 2 ](II) 
(gray). The platinum atoms are shown as white in the crystal structures, and light cyan in 
the averaged structures. 
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Figure A5. The averaged centroid structures of the platinum-DNA adduct and an 
undamaged B-DNA for: (a) 26.28DPZM-DNA; (b) 4.7DPZM-DNA; and (c) 3.7DPZM-DNA. 
[Pt2(DPZM)Cl 2 ](II) is shown as pink, bound-DNA as blue, and B-DNA as gray. 




(a) (b) (c) 



Figure A6. Overlap of the averaged structures for B-DNA (gray), [Pt2(DPZM)Cl2](II) and 
[Pt 2 (DPZMCH2)Cl 2 ](II) (pink), and: (a) 4.7DPZM-DNA (blue) and 4.7DPZMCH2-DNA 
(cyan); (b) 4.7DPZM-DNA (blue) and 6.28DPZM-DNA (royal); and (c) 4.7DPZMCH2-DNA 
(cyan) and 6.28DPZMCH2-DNA (green). 
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Figure A7. The average centroid structures of the 12 platinum-DNA adducts 
studied (pink for [Pt 2 (DTBPA)Cl 2 ](II) and [Pt 2 (TPXA)Cl 2 ](II); blue for intra-strand 
[Pt 2 (DTBPA)Cl 2 ](II)-bound DNA; navy for inter-strand [Pt 2 (DTBPA)Cl 2 ](II)-bound 
DNA; cyan for intra-strand [Pt 2 (TPXA)Cl 2 ](II)-bound DNA, green for inter-strand 
[Pt 2 (TPXA)Cl 2 ](II)-bound DNA; gray for B-DNA): (a) 26.28DTBPA-DNA; 
(b) 4.7DTBPA-DNA; (c) 3.7DTBPA-DNA; (d) 7.28DTBPA-DNA; (e) 6.28DTBPA-DNA; 
(f) 7.26DTBPA-DNA; (g) 26.28TPXA-DNA; (h) 4.7TPXA-DNA; (i) 3.7TPXA-DNA; 
(j) 7.28TPXA-DNA; (k) 6.28TPXA-DNA; (1) 7.26TPXA-DNA. 
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Figure A8. Selected frequency distributions of the representative DNA duplex base-pair 
helical parameters and base-pair step parameters for the central binding base-pairs of 
B-DNA (black), 26.28DTBPA-DNA (green), 3.7DTBPA-DNA (blue), 7.28DTBPA-DNA 
(dark cyan), 6.28DTBPA-DNA (magenta) and 7.26DTBPA-DNA (dark yellow); and 
26.28TPXA-DNA (red), 4.7TPXA-DNA (blue), 3.7TPXA-DNA (dark cyan), 7.28TPXA-DNA 
(magenta), 6.28TPXA-DNA (dark yellow) and 7.26TPXA-DNA (navy). 
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Figure A8. Cont. 
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Figure A9. Major groove widths and depths for the time-averaged structures of the DNA 
conformations of the other 1 1 platinum-DNA adducts. Major groove width (a) and depth 
(b) for [Pt 2 (DTBPA)Cl 2 ](II) are shown for: B-DNA (black squares); 26.28DTBPA-DNA 
(green circles); 3.7DTBPA-DNA (blue circles); 7.28DTBPA-DNA (dark cyan circles); 
6.28DTBPA-DNA (magenta circles); and 7.26DTBPA-DNA (dark yellow circles). Major 
groove width (c) and depth (d) for [Pt 2 (TPXA)Cl 2 ](II) are shown for: B-DNA (black 
squares); 26.28TPXA-DNA (red triangles); 4.7TPXA-DNA (blue triangles); 
3.7TPXA-DNA (dark cyan triangles); 7.28TPXA-DNA (magenta triangles); 
6.28TPXA-DNA (dark yellow triangles); and 7.26TPXA-DNA (navy triangles). 
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